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INTRODUCTION 

The  equations  for  the  stream  functions  of  a  steady,  two- 
dimensional  flow  of  a  compressible  fluid,  which  are  non¬ 
linear  when  considered  in  the  physical  plane,  become  linear 
when  transformed  to  the  hodograph  plane,  see  Ref.  11.  In 
this  latter  plane,  therefore,  the  principle  of  superposi¬ 
tion  can  be  applied.  This,  as  well  as  the  fact  that  the  ho¬ 
dograph  method  is  more  suitable  for  the  determination  of 
pressure  distribution  on  the  boundary  of  the  profile,  sug¬ 
gests  that  the  hodograph  method  be  used  in  developing  pro¬ 
cedures  for  the  computation  of  steady,  two-dimensional 
flows. 

The  procedures  for  the  computation  of  flow  patterns  by 
use  of  the  hodograph  method,  especially  the  transformation 
back  to  the  physical  plane,  require  a  very  considerable 
amount  of  numerical  work  even  in  the  case  of  an  incompres¬ 
sible  fluid.  Moreover,  the  problem  of  determining  the  com¬ 
plex  potential  around  a  profile  prescribed  in  the  physical 
plane  leads  to  very  complicated  non-linear  boundary  condi¬ 
tions  when  the  hodograph  method  is  used.  There  are  at  pres¬ 
ent  no  methods  available  for  solving  this  problem  directly. 
This  makes  it  desirable  to  prepare  tables  based  on  the 
principle  of  the  hodograph  method  but  yielding  the  image  of 
the  streamlines  directly  in  the  physical  plane.  Thus  the 
mathematical  simplicity  of  the  hodograph  method  could  be 
exploited,  and  at  the  same  time  the  bulk  of  the  numerical 
work  inherent  in  its  use  could  be  avoided.  In  the  present 
paper  a  method  for  constructing  such  tables  will  be  de¬ 
scribed. 

In  order  to  Illustrate  the  principle  along  which  such 
tables  may  be  prepared,  the  preparation  of  the  tables  for 

*  Research  paper  done  under  Navy  Contract  NOrd  8555-Task  F  at  Harvard 
University.  The  ideas  expressed  in  this  paper  represent  the  personal 
view  of  the  author  and  are  not  necessarily  those  of  the  Bureau  of  Ord¬ 
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an  Incompressible  fluid  will  first  be  explained.  By  making 
use  of  a  correspondence  principle  Introduced  in  Refs.  1,  3s 
7,  and  8,  between  incompressible  and  compressible  flow, 
this  method  will  be  extended  to  a  compressible  fluid.  Then 
two  alternative  methods  will  be  discussed  for  setting  up 
the  tables,  and  a  sample  will  be  presented,  prepared  ac¬ 
cording  to  one  of  these.  A  numerical  example  illustrating 
the  use  of  the  tables  will  then  be  given.  Finally,  the  con¬ 
ditions  under  which  a  given  hodograph  corresponds  to  a 
physically  possible  flow  pattern  will  be  discussed. 

The  author  wishes  to  thank  Dr.  Bernard  Epstein  for  val¬ 
uable  suggestions  and  Mr.  Maurice  Neuman  for  help  in  the 
preparation  of  the  present  paper. 

NOTATION 

a  -  local  velocity  of  sound,  dependent  upon  the  speed  of 
the  fluid 

ao  -  velocity  of  sound  at  the  stagnation  point 
k  -  ratio  of  isobarlc  to  Isovolumic  specific  heat 
p  -  pressure 

Po  -  pressure  at  the  stagnation  point 
q  *  -  speed 

x,y-  Cartesian  coordinates  in  the  physical  plane 
M  -  Mach  number  (M  =  q/a) 

0  -  angle  subtended  by  the  velocity  vector  with  the  x-axls 

p  -  density 

Po  -  density  at  the  stagnation  point  which  will  be  taken 
equal  to  1 

G  -  a  constant  in.  the  pressure  density  relation  (p  =  crp^) 

9  -  potential  function 

i|/  -  stream  function 

Whenever  a  function  is  expressed  in  terms  of  two  sets  of 
variables  obtained  from  each  other  by  point  transformations, 
the  same  symbol  will  be  retained  for  both  functional  forms. 
Thus  i|/(x,y)  =\j;(q,0)  when  x  =  x(  q,  0  )  and  y  ~  y  (  q^  0 )  • 

In  the  present  paper  flow  patterns  will  be  considered 
which  are  characterized  by  three  different  equations  of 
state: 

(1)  p  =  Po  =  const  =  1  for  an  incompressible  flow, 

(2)  p  =  (p/a)^/^  for  an  exact  compressible  flow, 

(3)  an  approximate  form  of  (2)  for  a  simplified  compres¬ 
sible  flow.  (See  the  second  section.) 
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Quantities  pertaining  to  incompressible  flow  Will  be  denot¬ 
ed  by  letters  capped  with  tlldas,  e.g.,  H,  ij/,  quanti¬ 

ties  pertaining  to  exact  compressible  flow  will  be  denoted 
by  plain  letters  (except  for  superscripts),  e.g.,  H,  X 
those  referring  to  simplified  compressible  flow  by  letters 
with  asterisks,  e  .g. ,  ^X*  etc. 

ON  TABLES  FOR  AN  INCOMPRESSIBLE  FLUID 

A  plane  is  introduced  whose  Cartesian  coordinates  are 
the  variables 


n  =/^Po((iq/q)  =  log  q  (l) 

and  8  (see  notation). 

The  complex  potential . for  a  large  class  of  symmetric  ob¬ 
stacles  in  these  variables  can  be  written  in  the  form 


00 

s(z)  =  ^  ~ 

V  =  o 


(2) 


with  z  -  ?]  +  19  and  rio  the  logarithm  of  the  speed  at  infin¬ 
ity,  In  the  case  of  symmetric  obstacles  the  b^  s  are  purely 
imaginary  (see  Ref.  10,  §  4). 

Remark :  The  representation  (2)  holds  only  in  a  circle 
with  origin  at  rj  =  no  and  radius,  the  distance  from?)  -  no  to 
the  nearest  singularity.  However,  by  employing  suitable 
summation  methods,  e.g.. 


g(z)  =  lim  X  [  bv/l'U+ )J  [z  -  (2a) 

a  representation  is  obtained  for  g(z)  in  the  largest  doubly 
covered  star  domain  in  which  the  function  is  regular.  (For 
details,  see  Ref.  12.)  With  the  b^s  purely  imaginary,  the 
stream  function  ij)  can  be  written 


Im[g(z)]  =  ^  lb  5('')  {fi,e; 


no 


V  =  o 


( n^  0;no)  =  Re[z- noJ''"^^^ 


(3) 
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The  streamlines  of  the  flow  patterns  corresponding  to  g(z) 
In  the  physical  plane  are  given  by 

CX)  “ 

X  =  ^  x('’>  (5),  e,Tio),  y  =  ^  Ib^Y^'’)  (f),e;r)o) 

M  =  o  v=  o  (4) 

where 

J},0;r)o)  =/(Vpoe^°)  i[-cos  -  sin  eii/^'')]dn 

+  [cos  -  sin  0ij)  d0}  (5a) 

“H  0 

ri,,0;r]  o)  =f  ( l/poe^°)  i  [  -  sin  0$  +  cos  05pi''^j  dr) 

®  n 

+  [sin  +  cos  d0  !  (5b) 

r\  0 

b^s  depend  on  the  boundary  curve  of  the  profile.  In  employ¬ 
ing  the  hodograph  method  In  the  incompressible  fluid  case, 
it  Is  convenient  to  prepare  a  set  of  tables  for 

$  ^''^  ( rji  0;  r)o)i  ( 7),  0;r)o),  ( r),  0;rio)  (6) 

for  a  number  of  values  of  the  parameter  r)o.  For,  by  using  a 
set  of  b^s  appropriate  to  the  obstacle,  $(?),0;j)o)  could  Im¬ 
mediately  be  determined  by  making  use  of  Eqs.  (5)  and  (6). 
To  find  the  Image  of  a  streamline 

$(?i^0;no)  =  s  =  const  (7) 

in  the  physical  plane,  a  sufficiently  dense  set  of  points 
(Jjg^eg)  =  1,2,...,  Is  first  found,  corresponding  to  the 
streamline  $(  Ji,0  jr]o )  =  s.  From  the  tables  values  are  read 
for 

=  x(''hfis.esJno),  (8) 

and  by  using  Eq.  (4)  the  Image  ij)(x,y)  =  s  of  ii;(J),0;7)o)  =  s 
is  found  In  the  physical  plane. 

The  essence  of  this  method  Is  the  separation  of  the 
quantities  x(v)  ,  y(v)  ,  iji(v)  that  are  independent  of  the 
profile  from  the  b^s  which  depend  on  it.  For  the  former,  a 
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set  of  tables  can  be  prepared  once  and  for  all.  The  prob¬ 
lem  of  determining  the  latter  recurs  with  each  profile. 

This  Is  a  relatively  simple  matter  In  the  case  of  incom¬ 
pressible  flow.  New  variables  rj,  6  are  Introduced  Into  the 
complex  potential  g(z)  and  an  expansion  is  made  about 
The  coefficients  of  the  expansion  are  the  desired*  j b  .  In 
the  compressible  case  the  problem  Is  much  more  Involved, 
but  In  most  cases  of  Interest  the  b^s  found  for  the  incom¬ 
pressible  case  may  be  used  In  the  compressible  case  as  a 
first  approximation  to  the  exact  b^s.  Methods  are  known  for 
Improving  on  this  approximation  but  this  matter  will  not  be 
discussed  in  the  present  paper. 

EXTENSION  OP  THE  METHODS  USED  IN  THE  CASE  OF  AN 
INCOMPRESSIBLE  FLUID  TO  THE  CASE  OF  A 
COMPRESSIBLE  FLUID 

The  complex  potential  for  the  incompressible  fluid  flow 
considered  in  the  fj^e  plane  has  (in  general)  singularities, 
such  as  branch  points,  poles,  etc.  The  procedures  developed 
for  the  treatment  of  this  case  make  use  of  various  proper¬ 
ties  of  analytic  functions  defined  on  a  Riemann  surface. 

The  recently  developed  theory  of  the  Integral  operators 
permits  the  extension  of  these  procedures  to  theory  func¬ 
tions  satisfying  general  partial  differential  equations,  es¬ 
pecially  those  satisfying  the  compressibility  equation.  (For 
details  see  Refs.  3,  5,  7,  8.)  In  particular,  a  corre¬ 

spondence  principle  can  be  set  up  between  the  flow  pattern 
of  an  incompressible  fluid  around  a  profile  C,  and  that^of 
a  compressible  fluid  around  a  profile  C,  approximating  C. 

It  is  easiest  to  display  this  correspondence  principle  if 
the  equations  for  the  compressible  fluid  case  are  trans¬ 
formed  to  the  ’’quasi  logarithmic"  plane  defined  by  the  Car¬ 
tesian  coordinates  (?i,0) 

ri(q)  =  n(q)  =/^p(<l)[d'3/q]  (9) 

with  the  constant  of  Integration  so  chosen  that  r|(a)  =  0**, 

*For  the  construction  of  flow  patterns  of  an  incompressible  fluid 
around  a  given  profile  by  the  use  of  the  hodograph  method,  see  Ref.  l4. 

**It  is  noted  that  in  considering  purely  subsonic  flows  in  previous 
papers  the  variable 

X  =  J  )dg,0 

was  used  and  for  purely  supersonic  flovs  the  variable  X=  1X,0.  In 
studying  transonic  flow  patterns  it  is  found  to  be  more  convenient  to 
use  the  variables  r^,0. 
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Thus  r|  <  0  and  rj  >  0  for  the  sub-  and  supersonic  ranges, 
respectively.  In  the  incompressible  case 

ri  (q)  =  r,  (q)  =  /'^  p  (dq/q)  =  (dq/q)  (9a) 

where  the  constant  of  Integration  Is  now  so  chosen  that 
r](q)  =  log  q.  The  transition  from  (J5,0)  to  (rj,0)  means  only 
a  distortion  in  the  n  direction.  In  this  coordinate  system 
the  compressibility  equation  becomes 

L(\}/)  =  +  't (rj ) =  0  (lO) 


^(n)  =  [1  -  M^(r,)j/p  =  (r)) 

M  Is  the  local  Mach  number,  M  =  q/a(q).  (See  Ref.  3,  Par,  7.) 
For  the  Incompressible  case,  M  =  0,  p  =  po  =  l  and  Eq.  (lO) 
reduces  to 


L(i|j)  =  +  ^^\jr/^0^  =  0 


(10a) 


From  the  expression  for  'i(r])  it  Is  seen  that  the  equation 
Is  elliptic  In  'the  subsonic  and  hyperbolic  in  the  supersonic 
case. 

By  assuming  an  adiabatic  pressure  density  relation 
p  =  crp^  and  employing  the  Bernoulli  equation,  the  following 
explicit  expression  for  a(q)  and  p(q)  Is  found: 

p(q)  =  Po  !l  -  ['(k  -  l)/2j  (q/a)^r/^'''^^ 

y  (11) 

a^(q)  =  ao^(q)  -  [(k-l)/2jq^ 

J 

(See  Ref.  1,  p.  8-13.)  Use  of  these  formulas  yields  after  a 
lengthy  computation  the  following  expression  for^(r])  (see 
Ref.  9^  appendix). 


^(r,)  =  2[2/(k+l)J^^-'^’/(k+^)j  -(2r]) 


-  [(3k+5)/(2k+2)][2/(k+l)]"V('"''‘\-27if  +  ...  i  (12) 

Eq.  (lO)  is  now  simplified  by  retaining  only  the  linear  term 
in  7]  of  Eq,  (12),  i.e., 

-€(0)  =  =  Cri,  C  =  -2[2/(k  +  l)]  (12a) 


Determination  of  Transonic  Flov  Patterns 


19 


and  writing 

L*(i|/ )  =  =  0  (lO*) 

Remark :  Tricomi  (see  Ref.  15)  was  the  first  to  study 
Eq.  (lO*)  extensively.  Different  authors  subsequently  con¬ 
sidered  it  in  connection  with  the  theory  of  compressible 
fluids.  However,  the  present  treatment  is  different  from 
the  treatments  given  previously.  The  "correspondence  prin¬ 
ciple"  is  introduced  and  solutions  of  this  equation  defined 
on  a  Rlemann  surface  and  possessing  singularities  are  con¬ 
sidered.  Eq.  (10*)  will  henceforth  be  referred  to  as  the 
"simplified  compressibility  equation," 

The  application  of  the  operator  method  to  the  solution 
of  Eq.  (10*)  will  now  be  described.  Every  analytic  function 
g  of  one  complex  variable  z  =  r]  +  16  can  be  written  in  the 
form 

g(z)  =  m[f(z)J  =  {  f[z(l-  t=)/2]/(l-  (lla) 

with 

f(z)  =  n[g(z)J  =  (2/7t)ig(0)  +  z  g'[2z(l-  t=)Jdti  (lib) 

where  g'(0  =  •  I't  has  been  shown  in  Ref.  9  that 

for  every  equation  of  type  (10*)  it  is  possible  to  define 
an  operator  p*  acting  on  g(z)  such  that 

U*(z)  =  P*[g(z)j  =  p*  {  m[f(z)]  1  =  p*[f(z)j  (13) 

is  a  complex  solution  of  Eq.  (10*).  The  following  are  the 
explicit  expressions  for  P*[f(z)J  and  p*[g(z)}; 

P*[f(z)]  =/^'-E*(ri,ejno)f[(z/2)(l-  )J  [( dt/  /l  -  t^  ]  (l4) 

P*[g(z)j  =  So(-Ti)"^/^  {  g(z)  + 

+  L  (l/2^^)[r(2n+  l)/r(n+  l)J  ^^(z)g^^^(z)  i  (15) 


'*”''Tliis  simplification  is  equivalent  to  replacing  the  exact  pressure 
density  relation  by  an  approximate  one.  It  is  possible  to  shov  that  the 
approximation  to  the  exact  equation  of  state  is  quite  good  in  a  certain 
range,  in  particular  around  the  transonic  line. 
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E*(r],0jj]o)  may  be  taken  as 

E*(r),e;no)  =  So(-ri)'^/^  «il/6,  5/6,  1/2, 

+  le  -  lno'\/2c{-^)^/^\  (16) 

where  ^(a^p^y^X)  is  the  hypergeometrlc  function.  So  and  c 
are  numerical  constants 

Sq  ^  2  (k  +  1)  (2-k) /(ek-e  ) 

0  =  (2VV5)[(k-l)/2](^‘'-®’/^®''-®^ 

The  polynomials  |i  (z)  are  defined  by  the  recursion  formu- 
la  and  Initial  conditions: 

+  2P*  =  0  *>1 

(2n  +  l)l)!^^“'^^Mz) +2?^!"'  +2P*?”(z)=0  >•  (18) 

t3(")(0)  =0  J 

with  F*  =  (5/36)  (l/4o^)  (-rj)-^ 

The  functions  may  be  written  as 


^g(zn)  az=  [(-l)7(n-l):]  r 


^(z-t)^  g(t)dt 


The  method  of  Integral  operators  permits  setting  up  a 
correspondence  principle  between  the  flow  patterns  of  an  In¬ 
compressible  and  a  compressible  fluid.  Both  fluids  are  con¬ 
sidered  In  the  quasi  logarithmic  plane  (which  In  the  case 
of  an  Incompressible  fluid  reduces  to  the  logarithmic  plane). 
The  correspondence  principle  consists  In  associating  with 
the  flow  pattern  ij-  =  Im[g]of  an  Incompressible  fluid  a  flow 
pattern  \{i  =  Im  [p*(g)]  of  a  compressible  fluid.  By  intro¬ 
ducing  expansion  (2)  for  g(z)  expressions  are  obtained  for 
a  compressible  fluid  that  are  the  analogues  of  Eqs.  (3)> 

(4),  and  (5)  for  an  Incompressible  fluid.  Thus 

p*[g{z)]  =  Z  7P*[(z-rio)''"'-/^J  (20) 

where 

P*[(z  -  no)'"^^'^^^  =  p*(z),  p*(z);)lo][p*(z)  -  11  (21) 
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and 

p*  (z)  =  p*[ri  *t  10j  =  -c  +  10  =C  (22) 

b('')  [p*(z),p*(z);r,o]  =  B<'')  [!:,C;no] 

=  So(-»,)-"/‘‘«[l/6,5/6,v+l/2,(C-r,o)/(-2c)(-r,)Va]  (23) 

It  is  seen  then  that  in  the  case  of  the  stream  function  the 
transition  from  the  Incompresslhle  to  the  compressible  case 
is  made  by  multiplying  each  term  of  its  power  series  expan¬ 
sion  by  a  "compressibility  factor"  B  (^) ( Cj  ?lo)  and  by  re¬ 
placing  the  variable  z  =  r]  +  10  by  the  variable 
2;  =  _c(-rj)V®  +  10.  This  latter  change  implies  a  distor¬ 
tion  in  the  rj  direction  mentioned  above. 

Corresponding  to  expression  (3)  of  the  Incompressible 
case,  the  stream  function  for  a  symmetric  profile  can  be 
written  in  the  form 

=  X  1 

where 

(25) 

The  implicit  representation  of  the  streamline  h*  (^1^  ’lo)  =  s 
in  the  physical  plane  is  given  by 

X  =  r  X<'’hn»05ilo),  y  =  L  ib^  y^'’^  (r),e;7io)  (26) 

with 

X*’’  =;  )  ( [_{1  .  cosB-<,y  slneldi, 

+  AT'Iin 003 e  -  ifg’*  Slnelde]  (27a) 

y'*’  =/  (l/p)(l/e^  sine-,''')  coaeldp 

+  sine+  ^y’oosejdet 


(27b) 
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In  order  to  evaluate  the  "compressibility  factors" 

B  (^)  ( r],0;rjo) ,  suitable  power  series  representations  must  be 
found  for  the  hypergeometrlc  function  appearing  in  Eq.  (23). 
It  is  well  known  that  for  |  x|  <1,  repre¬ 
sented  by  a  hypergeometrlc  series  P(  a,  and  for 

I  X  I  >  0  a  combination  of  two  hypergeometrlc  series  with 
arguments  and  parameters  different  from  those  appearing  in 
the  hypergeometrlc  function  must  be  used  (see  Ref.  l6,  §  IV). 
In  the  case  of  the  hypergeometrlc  function  appearing  in  Eq, 
(23),  the  boundary  separating  the  domains  of  convergence, 

-  ,50  +  ieJ/[-2c{-r|)^/^]  I  =  1 


e 


Fig.  1,  Domains  of  Validity  of  the  Power  Series  Representations 
(28a)  and  (281))  for  the  Hypergeometrlc  Function 

-c(-r,)^/®  -  rio  +  ie\ 

(A)  3c®r®  +  0®  +  2c(-i))®/^  n  O  +  =  0  <  0 

(B)  3o^r)®  -  0®  -  2c(-r))®/^  0  _  ^  q  y,  >  q 

In  Si  =  si  +  Si  and  in  Sg  =  S2  +  S2,  the  representations 

^[6’6’  ^2’  .2c(-nf/^  )  =^[6’6’  ^2'  .2o(-r,f/^  I 

(28a) 


V  + 
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and 


,o(-o)^/^  -  r,0  +  le 


■’6’  ^  2’ 

fl2 

-2c(-r))^/^-  r)  0  +  10^ 

)  ' 

D  p 

V 

5/6 

F 

5 

-2o(-r,)®/®-  no  +ieJ 

6’3 

-2c(-ri) 


3/2 


-c  ( -  Tj)  r]  o  +  10 


(28b) 


are  used  respectively.  The  A^s and  BJs  denote  the  quantities 
\  =  r(-2/3)r(v  +  i/2)/r(-i/3  -  v  )(5/6) 

=r(v  +  i/2)r(2/3)/r(i/3  -*v)r(i/6) 


TABLES  FOR  A  COMPRESSIBLE  FLUID 

It  has  been  shown  in  the  previous  section  how  it  is  pos¬ 
sible  to  extend  the  methods  used  in  the  case  of  an  incom¬ 
pressible  fluid  to  the  case  of  a  compressible  one.  The  an¬ 
alogues  of  Eqs.  (2),  (3),  (4),  and  (5)  of  the  incompressi¬ 
ble  flow  were  obtained,  namedly  Eqs.  (20),  (25)^  (26),  and 
(27).  These  permit  the  setting  up  of  tables  for  a  compres¬ 
sible  fluid  of  the  type  that  was  outlined  in  the  first  sec¬ 
tion  for  an  Incompressible  fluid. 

In  this  section  two  types  of  tables  will  be  indicated. 
Type  I  depends  upon  a  parameter  r]  o  and  has  to  be  prepared 
for  a  number  of  values  of  rj  sayr]o,v(’^  =1,  , , ,  n) ,  In 
this  case  one  must  resort  to  interpolation  methods  in 
order  to  obtain  streamlines  for  Intermediary  values  of 
ilo  :  ^  o>v  <  ’lo  <  ilo^v+i«  Tables  of  type  II  are  Independent 
of  any  parameter  but  their  application  requires  more  numer¬ 
ical  work  when  used  for  a  value  of  ^  o,  for  which  tables  of 
type  I  are  available. 

Type  I,  In  the  present  section  a  table  of  (r] , 6  o) 

for  a  specific  value  of  =  .-#02  will  be  given.  A  table 

is  also  included  for  . 

j-.02)  =  Imip*[  (z  -  .02)''~^  i 
in  view  of  possible  later  applications,  and  a  conversion 
table  from  the  variable  r)  to  the  local  Mach  number  M.  An  il¬ 
lustration  of  the  use  of  these  tables  will  be  given  in  the 
next  section. 
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TABLE  I.  VALUES  OF  THE  FUNCTIONS 


,(v) 


(Jl.  e) 


0 

V  =  0 

V  =  1 

V  =  2 

V  =  3 

V  =  4 

V  =  5 

-  .050 

0.5 

1.876 

.015 

-.869 

-.043 

.272 

.025 

1.0 

1.003 

.012 

-2.049 

-.105 

2.457 

.209 

1-5 

.802 

.008 

-2.985 

-.161 

8.419 

.757 

-.045 

0.5 

1.859 

.015 

-  .854 

-.043 

.265 

.024 

1.0 

1.010 

.010 

-2.026 

-.101 

2.4o4 

.214 

1.5 

.798 

.008 

-2.951 

-.164 

8.577 

.772 

-.040 

0.5 

1.842 

.016 

-.838 

-.042 

.258 

.025 

1.0 

1.017 

.010 

-2.002 

-.100 

2.407 

.217 

1.5 

794 

.008 

-2.920 

-.165 

8.696 

.784 

-.035 

0.5 

1.825 

.016 

-.824 

-.04l 

.252 

.023 

1.0 

1.024 

.010 

-1.979 

-.099 

2.562 

.213 

1.5 

.790 

.008 

-2.888 

-.167 

8.832 

.797 

-.030 

0.5 

1.809 

.016' 

-.809 

-.040 

.245 

.022 

1.0 

1.030 

.010 

-1.954 

-.098 

2.315 

.208 

1.5 

.785 

.008 

-2.857 

-.168 

8.929 

.808 

-.023 

0.5 

1.790 

.016 

-.793 

0 

0 

.257 

.021 

1.0 

1.037 

.010 

-1.929 

-.096 

2.267 

.204 

1.5 

.781 

.008 

-2.825 

-.170 

9.115 

.821 

-.020 

0.5 

1.771 

.016 

-.778 

-.039 

ro 

0 

.021 

1.0 

1.044 

.010 

-1.905 

-.095 

2.221 

.200 

1.5 

.777 

.008 

-2.794 

-.171 

9.252 

.838 

-0.015 

0 

9.669 

0.157 

0.002 

_ * 

— 

— 

0.5 

1.754 

0.017 

-0.763 

-0.038 

0.223 

0.020 

1.0 

1.051 

0.011 

-1.88 

-0.094 

2.173 

0.196 

1.5 

0.805 

0.007 

-3.012 

-0.156 

8.274 

0.734 

-0.010 

0 

8.667 

0.179 

0.003 

— 

— 

— 

0.5 

1.737 

0.017 

-0.748 

-0.037 

0.215 

0.019 

1.0 

1.058 

0.011 

-1.857 

-0.092 

2.126 

0.192 

1.5 

0.810 

0.007 

-3.045 

-0.157 

8.139 

0.721 

-0.005 

0 

7.927 

0.198 

0.004 

— 

— 

— 

0.5 

1.715 

0.017 

-0.732 

-0.037 

0.208 

0.019 

1.0 

1.063 

0.011 

-1.830 

-0.092 

2.076 

0.187 

1.5 

0.813 

0.0C8 

-3.071 

-0.155 

7.985 

0.715 

0.005 

0 

13.222 

0.472 

0.012 

— 

— 

— 

0.5 

0.897 

1.125 

0.318 

-0.359 

0.072 

0.102 

1.0 

0.551 

1.4o6 

-0.841 

-1.761 

0.872 

1.963 

1.5 

0.416 

1.6o4 

-1.492 

-4.454 

3.350 

11.146 

0.010 

0 

11.839 

0.516 

0.013 

— 

— 

— 

0.5 

0.906 

1.138 

-0.310 

-0.369 

0 .068 

0.106 

1.0 

0.555 

1.417 

-0.837 

-1.785 

0 .847 

2.014 

1.5 

0.408 

1.615 

-1.474 

-4.522 

3.479 

11.384 

The  dash  indicates  that  the  value  is  smaller  than  ,0005 
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TABLE  I. 


VALUES  OF  THE  FUNCTIONS 


(J1  ^  6) 


(Cont. ) 


n 

8 

V  =  0 

V  =  1 

V  =  2 

V  =  5 

V  =  4 

V  -  5 

0.015 

0 

10.357 

0.556 

0,014 

* 

0.5 

0.914 

1.149 

-0.302 

-0.379 

0.064 

0.110 

1.0 

0.557 

1.423 

-0.824 

-1.811 

0.819 

2.058 

1.5 

0.422 

1.622 

-1.458 

-4.585 

3.4o6 

11.600 

0.020 

0 

10.296 

0.599 

0.017 

— 

— 

— 

0.5 

0.922 

1.161 

-0.294 

-0.388 

0.059 

0.113 

1.0 

0.561 

1.434 

-0.813 

-1.844 

0.794 

2.108 

1.5 

0.423 

1.63 

-1.435 

-4.644 

3.333 

11.818 

♦  025 

0.5 

.931 

1.173 

-.286 

-.598 

.055 

.117 

1.0 

.565 

1.445 

-0.802 

-1.876 

.769 

2.159 

1.5 

.426 

1.639 

-1.426 

-4.704 

5.260 

12.032 

.030 

0.5 

.944 

1.185 

".279 

-.407 

.051 

.121 

1.0 

.569 

1.454 

-.790 

-1.908 

.746 

2.208 

1.5 

.428 

1.650 

-1.411 

-4.766 

3.187 

12.255 

.035 

0.5 

1.197 

-.270 

-.417 

.046 

.125 

1,0 

.572 

1.464 

-.777 

-1.959 

,720 

2.258 

1.5 

.430 

1.656 

-1.395 

-4.827 

3.113 

12.474 

.040 

0.5 

.955 

1.209 

-.263 

-.427 

.043 

.129 

1.0 

.575 

1.473 

-.765 

-1.969 

.695 

2.508 

1.5 

.432 

1.664 

-1.579 

-4.896 

5.041 

12.693 

.045 

0.5 

.964 

1.220 

-.255 

-.437 

.038 

.155 

1.0 

.578 

1.482 

-.753 

-2.000 

.670 

2.359 

1.5 

.454 

1.66 

-1.363 

-4.953 

2.968 

12.918 

.050 

0.5 

.973 

1.251 

-.247 

-.446 

.055 

.137 

1.0 

.582 

1.491 

-.741 

-2.051 

.645 

2.410 

1.5 

.i^■33 

1.671 

-1.547 

-5.009 

2.896 

15.142 

*  The  dash  indicates  that  the  value  is  smaller  than  .0005 
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(  V  )  ,  , 

T.'ffiLE  II,  THE  VALUES  OE  FUNCTIONS  ;  0 ) 


0 

V  =  0 

V  =1 

V  =  2 

V  =  3 

V  =  4 

V  =  5 

-0.050 

0.5 

.017 

.668 

.020 

-.234 

-.017 

.069 

1.0 

.012 

1.688 

,051 

-2.299 

-.162 

2.686 

1.5 

.007 

1.935 

.053 

-5.697 

-.281 

14.761 

-  .045 

0.5 

.017 

.644 

.019 

-.223 

-.016 

.065 

1.0 

.015 

1.706 

.051 

-2.260 

-.153 

2.619 

1.5 

.007 

1.926 

.053 

-5.628 

-.292 

14.517 

-.040 

0.5 

.017 

.619 

.019 

-.212 

-.015 

.062 

1.0 

.013 

1.694 

.051 

-2.224 

-.156 

2.564 

1.5 

.007 

1.916 

.053 

-5.558 

-.350 

14.268 

-.055 

0.5 

.017 

.593 

.018 

-.200 

-.014 

.058 

1.0 

.013 

1.684 

.051 

-2.190 

-.154 

2.510 

1.5 

.007 

1.907 

.054 

-5.488 

-.385 

14.019 

-.030 

0.5 

.016 

.565 

.017 

-.188 

-.013 

.054 

1.0 

.013 

1.673 

.050 

-2.154 

-.151 

2.453 

1.5 

.007 

1.901 

.054 

-5.419 

-.379 

13.769 

-  .025 

.  0-5 

.016 

.535 

.016 

-.180 

-.012 

.050 

1.0 

.015 

1.662 

.050 

-2.117 

i 

00 

2.395 

1.5 

.007 

1.888 

.055 

-5,346 

-.375 

15.515 

-  .020 

0.5 

.016 

.500 

.015 

~,l6l 

-.011 

.046 

1.0 

.014 

1.649 

.049 

-2.081 

-.145 

2.285 

1.5 

.007 

1.879 

.054 

-5.277 

-.369 

13.263 

-0.015 

0 

9.51T 

0.478 

0.013 

_ * 

— 

— 

0.5 

0.042 

1.317 

0.039 

-0.417 

-0.029 

0.117 

1.0 

0.013 

1.641 

0.047 

-2.045 

-0.143 

2.281 

1.5 

0.007 

1.840 

0.052 

-5.203 

-0.357 

13.000 

-0.010 

0 

10.416 

0.446 

0.011 

— 

— 

— 

1.5 

0.043 

1.303 

0.039 

-0.406 

-0,028 

0.112 

1.0 

0.013 

1.629 

0.049 

-2.009 

-0.141 

2.222 

1-5 

0.007 

1.860 

0.054 

.5.134 

-0.358 

12.732 

-0.005 

0 

11.478 

0.411 

0.009 

_ 

_ 

— 

0.5 

0.044 

1.289 

0.038 

-0.394 

-0.027 

0.107 

1.0 

0.014 

1.616 

0.048 

-1.970 

-0.138 

2.161 

1.5 

0.007 

1.847 

0.054 

-5.050 

-0.352 

12.470 

0.005 

0 

-0.002 

0.002 

— 

— 

— 

— 

0.5 

-1.432 

0 .616 

0.652 

-0.154 

-0.193 

0.0352 

1.0 

-0.914 

0.788 

1.601 

-0.868 

-1.866 

0.862 

1.5 

-0.700 

0.908 

2.739 

-2.325 

-7.086 

5.362 

0.010 

0 

-0.001 

— 

— 

— 

- : 

— 

0.5 

-1.4l8 

0.611 

0.667 

-0.148 

-0.199 

0.031 

1.0 

-0.909 

0.784 

1.631 

-0.852 

-1.908 

0.830 

1.5 

-0.687 

0.904 

2.773 

-2.288 

-7.218 

5.227 

The  dash  Indicates  that  the  value  is  smaller  than  .005, 
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TABLE  II.  THE  VALUES  OF  FUUCTIOUS  6)  (Cont. ) 


u 

0 

V  =  0 

V  =  1 

V  =  2 

V  =  5 

V  =  4 

V  =  5 

0.015 

0 

-0.002 

_ _ * 

0.5 

-l.ii-05 

0.6o4 

0.679 

-0.142 

-0.205 

0.027 

1.0 

-0.901 

0.777 

1.649 

-0.851 

-1.9i^5 

0.795 

1.5 

-0.696 

0.900 

2.800 

-2.252 

-7.577 

5.091 

0.020 

0 

-0.285 

0.007 

-0.001 

— 

— 

— 

0.5 

-1.586 

0.597 

0.692 

-0.155 

-0.211 

0.024 

1.0 

-0.895 

0.775 

1.672 

-0.8l4 

-1.985 

0.765 

1.5 

-0.691 

0.895 

2 .8220 

-2.215 

-7.456 

4.952 

.025 

0.5 

“I.57I 

.591 

.705 

-.129 

-.218 

.022 

1.0 

-.890 

.769 

1.695 

-.796 

-2.025 

.755 

1.5 

-.687 

.891 

2.854 

-2.179 

-7.574 

4.815 

.050 

0.5 

-1.550 

.585 

.719 

-.125 

-.224 

.019 

1.0 

-.88ii. 

.768 

1.716 

-.777 

-2.067 

.702 

1.5 

-.684 

.889 

2.882 

-2.145 

-7.694 

4.678 

.055 

0.5 

-1.540 

.578 

.751 

-.117 

-.250 

.017 

1.0 

-.877 

.758 

1.756 

-.758 

-2.108 

.670 

1.5 

-.680 

.882 

2.909 

-2.106 

-7.8i4 

4.540 

,0h0 

0.5 

-1.525 

.571 

.744 

-.111 

-.257 

.014 

1.0 

-.871 

.755 

1.757 

-.759 

-2.149 

.658 

1.5 

-.676 

.877 

2.956 

-2.067 

-7.955 

4.4o6 

.oil-5 

0.5 

-1.509 

.565 

.757 

-.105 

-.245 

.011 

1.0 

-.864 

.7i(-8 

1.778 

-.721 

-2,190 

.606 

1.5 

-.675 

.877 

2.965 

-2.018 

-8.054 

4.271 

.050 

0.5 

-1.295 

.558 

.771 

.099 

-.250 

.009 

1.0 

-.858 

.743 

1.798 

-.702 

-2.251 

.575 

1.5 

-.678 

.871 

2.990 

-1.997 

-8.175 

4.155 

*  The  dash  indicates  that  the  value  is  smaller  than  .005. 
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TABLE  III 


'"•onverslon  of  M  to  r] 


M 

n 

M 

.60 

-.1977 

.96 

-.0165 

.62 

-.1857 

.98 

-.0085 

.64 

-.1750 

1.00 

0 

.66 

-.1616 

1.01 

.0044 

.68 

-.15^1 

1.02 

.0079 

.70 

-.1407 

1.05 

.0120 

.72 

-.1507 

1.04 

.0160 

.74 

-.1195 

1.05 

.0197 

.76 

- . 1087 

1.06 

.0234 

.78 

-.0978 

1.07 

.0274 

.80 

-.0887 

1.08 

.0312 

.82 

-.0787 

1.09 

.0350 

.84 

-.0692 

1.10 

.0387 

.86 

-.0603 

1.10 

.0387 

.88 

-.0513 

.90 

-.0421 

.92 

-.0336 

.94 

-.0245 

Type  II.  In  preparing  this  type  of  table  the  values  of 
>((1/6,  5/6,  V  +  1/2,  a  +  Ip)  are  computed  once  and  for  all, 
as  well  as  Its  derivative  with  respect  to  the  argument  for 
a  number  of  values  of  a  and  p  where 

a  =  1/2  -  no/[-2c(-!i)^/®]  f!  =  0/[-2c(-n)^/®] 

By  the  use  of  tables  of  fractional  powers  and  trigonometric 
functions  and  the  tables  outlined  above,  the  computation 
can  now  be  made, 

=  Ee[So(-7i)]^/^[-c(-  rj)V^-  210  +  10]“^/^  X(1/6,5/6,v+ 1/2,  a  +  iP) 

and  (r],0;T]o),  ’^no)  * 

The  stream  functions  in  the  hodograph  plane  are  obtained 
by  forming  i}/(r],0;r]o)=  Zlby\(/(^)(r],0;r]o).  The  Inte- 

of  Eqs.  (27a)  and  (27b)  are  now  evaluated  by  using  the  ta¬ 
bles  for  the  derivative  of  X(l/6,  5/y>  v  +  1/2,  a  +  ip). 
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The  transition  to  the  physical  plane  Is  carried  out  by  nu¬ 
merical  Integration  of  Eqs,  (27a)  and  (27b), 

Note,  Since  the  integrands  of  Eq,  (17)  are  perfect  dif¬ 
ferentials,  the  value  of  the  integral  is  Independent  of  the 
path  and  it  may  conveniently  be  taken  along  a  streamline. 

Tables  of  Type  II  will  not  be  given  in  the  present  paper. 

ILLUSTRATION  OE  THE  USE  OP  THE  TABLES 

As  in  a  previous  paper^°  the  complex  potential  of  an  in¬ 
compressible  fluid  flow  around  an  elliptic  profile  is  chos¬ 
en  as  the  analytic  function  g(z).  As  is  well  known,  the 
complex  potential  of  an  incompressible  flow  around  a  circle 
is  given  in  the  physical  plane  (t  plane)  by  the  relation 

w  =  9  +  i^  =  qo(t  +  R/t^)  (29) 

where  qo  is  the  dimensionless  speed  at  infinity.  The  func¬ 
tion 

i  =  (1/2) (t  +  i/t)  (30) 


maps  the  exterior  of  the  circle  onto  the  exterior  of  the  el¬ 
lipse  and  therefore  the  function  g(0  =  w[t(^)]  will  repre¬ 
sent  the  complex  potential  of  an  incompressible  flow  around 
an  ellipse  in  the  physical  plane.  This  potential  when  trans¬ 
formed  to  the  hodograph  plane  will  furnish,  upon  expansion 
about  a  set  of  j  b^J  for  the  incompre'ssible  case,  which, 
as  a  first  approximation  will  also  be  used  for  the  compres¬ 
sible  flow.  (See  the  first  section.) 

This  transformation  is  carried  out  as  follows;  The  con¬ 
jugate  to  the  velocity  vector  q  is  given  by 

9  =  dw/d6=  (dw/dt)/(dg/at)  =  [(1/2)90(1 (1/2)90(1  -  1/t^)]  (31) 

or 


=  (q  -  qoR^)/(q  -  qo)  (32) 

Where  qo  is  the  velocity  at  infinity. 

Therefore, 

vig[t(9)]l,=  (l/2)9oU(q-qoR^)/Ci-9o)f/^  +  R®[(i-<lo)/(q-9oR^)]"^^! 


(33) 
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Since  the  complex  potential  In  the  logarithmic  plane  Is 

needed,  z  Is  set  =  rj  -  16  with  ^  =  J^dq/q  =  log  q  and 
T]o  =  log  qo.  Thus  Is  obtained 

q(z)  =  (l/2)Re(e'’°)  i[(l-R‘=e"  ’'°)/(l-e^ 

+  [(l-e'”'°)/(l-R-"e''''°)]^/"  i  (34) 

Which  can  be  developed  In  a  power  series  of  the  form 


g(z)  =  ^  bv(z  -  (34a) 

V  =  o 

As  an  example,  qo  is  taken  equal  to  O.98  (or  r)o  =  -0.02) 
and  =  1.2,  which  yields  a  thickness  ratio  for  the  ellip¬ 
tic  profile  of  11:1,  With  these  numerical  values  for  u  and 
R  the  set  of  (b^  !  listed  In  Table  IV  Is  obtained. 

TABLE  IV 

The  Series  Coefficients  of  Function  g 
V  0  12  3  ^  5 


b  .240101  -2.100861  -4.859411  -8.170151  -82.869231  -389.551 


It  Is  observed  that  the  radius  of  convergence  of  this 
series  Is  \  r]o  \  =  log  =  -18  about  the  center  ^;  =  r\o. 
As  was  mentioned  before,  however,  by  employing  suitable 
summation  methods,  the  series  may  be  used  outside  Its  do¬ 
main  of  convergence  (see  the  first  section). 

Using  these  values  for  {by!  and  the  table  of 
yields  the  expression 


>f  =  X  lb  -0.02) 

which  Is  equated  to  zero.  If  this  is  solved  by  a  graphical 
method,  a  number  of  points  are  obtained  on  the  streamline 

i(/(ri,e; -0.02)  =  0  (35) 


(  see  Table  V) . 
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TABLE  V 


Points  on  Streamline  \i/ (t)  ,  0;- .02 )  =  0 


r] 

.052 

.029 

.025 

.014 

0  - . 0l4 

0 

0 

.08 

12k 

.151 

170  .190 

The  next  task  Is  to  obtain  the  Image  of  this  streamline 
In  the  physical  plane.  To  do  this  the  values  of  and 

{Y  !  v=  0,1,  ...5,  are  found  corresponding  to  the  set 
(r|,e).  (See  Tables  VI  and  VII.) 

TABLE  VI 

( v) 

Values  of  X  (r]  ,0;- 0.02)  for  Points  on  Streamline  \|/ (t],  9;-0.02)  =  0 


n 

0 

,(0) 

z(5) 

z(5) 

.052 

0 

.26080 

.52965 

.02475 

.00054 

.00001 

0 

.029 

.08 

-7.L4056 

-2.19674 

2,10610 

.09445 

-.04775 

- .00019 

.025 

.124 

4.18929 

.92197 

.45728 

.15956 

-.00951 

.00272 

.014 

.151 

5.57784 

5.41555 

1.45748 

.59867 

-.15616 

.00155 

0 

.170 

.50167 

-.04515 

.02595 

.00582 

-.00119 

-.00052 

-.014 

.190 

8.00948 

-1.78920 

-1.86985 

-.44607 

.15081 

.05758 

TABLE  VII 

Values  of 

(r],e; 

-02)  for  : 

Points  on 

Streamline  \j/ (r(,0  j  - 

b 

ro 

II 

0 

n 

0 

y(0) 

,(1) 

y(2) 

^(3) 

y(5) 

.052 

0 

0 

0 

0 

0 

0 

0 

.029 

.08 

.46715 

- .17611 

.16885 

,00758 

-.00585 

-.00001 

.025 

.124 

.52218 

.11492 

.05700 

.01757 

-.00119 

.00054 

.014 

.151 

.54465 

.51988 

.21895 

.09112 

-.02577 

.00025 

0 

.170 

.05175 

-.00740 

.00445 

.00100 

-.00020 

-.00005 

-.014 

.190 

1.02992 

-.22806 

- .25916 

-.05751 

.01908 

.00481 
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By  forming  the  expressions 

X  =  2:  ibv  e;-0.02) 

y  =  X  Ibv  Y  e;-0.02) 

a  set  of  points  is  obtained  lying  on  the  streamline  ^  =  0 
in  the  physical  plane  (see  Table  VIII). 

TABLE  VIII 


Points  on  the  Streamline \j/  (ri,0;  -.02)  =  0  in  the  Physical  Plane 


.052 

.029 

.025 

.014 

0 

-.014 

0 

0 

.08 

.124 

.151 

.170 

.190 

M 

1.080 

1.072 

1.057 

1.035 

1 

.969 

X 

-2.71 

-4.68 

-10.71 

-21.68 

-21.78 

-25.25 

y 

4.204 

5.454 

1.594 

1.57^ 

1.084 

Fig,  2  represents  the  streamline  of  the  incompressible 
flow  ^(x,y)  =  0,  coinciding  with  the  boundary  curve  of  the 


Fig.  2.  Parts  of  the  Streamlines  =  ^  about  Oval¬ 

shaped  Obstacles  C  and  C  respectively. 
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obstacle  together  with  the  streamline  =  0  of  the 

compressible  flow.  In  setting  up  the  diagram,  the  two 
streamlines  were  drawn  in  different  scales  in  order  to  make 
the  terminal  points  of  the  two  ovals  coincide.  The  neigh¬ 
borhood  of  the  terminals  of  the  ovals  in  the  physical  plane 
corresponds  to  the  neighborhood  of  low  Mach  numbers  in  the 
hodograph  plane.  The  streamline  In  these  neighborhoods  In 
the  former  plane  Is  obtained  by  expanding  and 

Y  ;J|o)  in  terms  of  the  Mach  number  as  a  parameter  of 

smallness: 

x^'')(n,e;no)  =  x('’)(n,ejr)o)  +  m  Xi^'’\5i, ejrjo)  +  ... 


Y  +  m  +  ... 

Where  X(7i,0;rio)^  Y(r],0;r]o)  are  defined  by  Eqs.  (5a)  and 
{5b).  Only  the  part  of  the  streamline  located  in  the  upper 
half  of  the  physical  plane  has  been  plotted  In  the  diagram. 
The  part  In  the  lower  half  is  symmetric  to  It  with  respect 
to  the  x-axls. 

CONDITIONS  UNDER  WHICH  A  GIVEN  HODOGRAPH  CORRESPONDS  TO 
A  PHYSICALLY  POSSIBLE  FLOW  PATTERN 

In  the  considerations  of  the  previous  sections  the  ques¬ 
tion  was  completely  disregarded  as  to  whether  the  flow^  pat¬ 
tern  obtained  by  this  method  Is  of  physical  significance. 
Obviously,  If  the  function  H'(x,y)  Is  to  yield  a  flow  pat¬ 
tern  around  a  closed  profile  in  the  physical  plane.  It  is 
necessary  that 

(a)  the  streamline  'p(x,y)  =  0  (or  a  part  of  this  stream¬ 
line,  in  the  present  case  the  image  of  the  lines  ABA’B', 
see  Fig.  5)  yield  a  closed  curve  In  the  physical  plane; 

(b)  that  every  point  of  the  physical  plane  oe  covered  at 
most  once  (that  is,  that  the  mapping  unto  the  x,y-plane  be 
schllcht*). 

The  two  conditions  for  the  subsonic  case  will  now  be  con¬ 
sidered: 

*It  is  remarked,  however,  that  if  the  whole  image  in  the  x,y-plane 
is  not  schlicht,  the  flow  pattern  obtained  may  nevertheless  he  of  phys¬ 
ical  interest  since  in  some  instances  it  is  required  to  determine  flow 
patterns  covering  only  part  of  the  x,y-plane  or  only  a  part  of  the  flow 
pattern  obtained  is  used. 


Bergman 


54 


h 

Fig,  3,  The  Streamline  i|»  =  0  in  the  Hodograph  and  Physical  Planes, 

(a»)  The  necessary  and  sufficient  conditions  that  the 
boundary  AB  A*B»  when  mapped  by  Eq.  (26)  unto  the  physical 
plane  yield  a  closed  curve  In  the  latter  plane  have  been  es¬ 
tablished  in  a  previous  paper  for  flow  patterns  that  are'  not 
necessarily  symmetric.  (See  Ref.  7«) 

(b»)  The  determinant  of  the  mapping  from  the  hodograph  to 
the  physical  plane  ^ (x,y )/^ (u, v)  is  always  positive.  It  is 
therefore  sufficient  and  necessary,  in  order  that  the. map¬ 
ping  be  schllcht,  that  the  boundary  curve  AiBi  AiBi  be 
schllcht,  and  that  the  mapping  satisfy  certain  regularity 
conditions  at  infinity,  which  will  be  stated  presently. 

In  the  case  under  consideration  the  schllchtness  of  the 
boundary  curve  AiBiAIbI  follows  from  the  symmetry  of  the 
image  in  the  hodograph  plane  and  the  fact  that  the  hodograph 
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lies  in  the  half  plane  u  >  0.  Indeed,  if  the  point  P  moves 
along  a  line  AMP,  then  the  velocity  component'  u  is  always 
positive,  that  is,  the  x  coordinate  always  Increases.  The 
segment  BNA’  cannot  intersect  Itself  or  the  arc  ANB,  since 
the  coordinate  still  Increases  if  the  point  P  moves  along 
this  curve.  Prom  the  symmetry  properties  it  follows  that 
this  arc  lies  in  the  upper  half,  and  that  the  image  of 
A’N’B’A*,  in  the  lower  half  of  the  xy-  plane.  The  whole 
boundary  curve  is  then  schlicht.  The  mapping  must  satisfy 
the  following  conditions  at  Infinity:  As  a  doubly  covered 
circle  about  the  point  0  the  image  of  infinity  of  the  hodo- 
graph  plane  shrinks  to  a  point,  its  image  in  the  Inverted 
physical  plane  [x'  +  iy '  =  l/(x  +  iy) ]  must  be  schlicht 
and  also  shrink  to  a  point. 

In  generalizing  these  results  to  transonic  flow  patterns, 
one  must  be  on  guard  against  the  following  difficulties: 

(a)  The  functional  determinant  ^(u, v)/^(x,y)  is  not  nec¬ 
essarily  positive, 

(b)  The  mixed  character  of  the  equation  of  flow  may  give 
rise  to  hyperbolic  singularities  which  do  not  appear  in  the 
subsonic  case. 

The  general  character  of  mappings  of  this  type  has  been 
extensively  investigated  by  Morse  and  others  (see  Ref.  1^) . 
This  problem  will  not  be  discussed  in  the  present  paper. 
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